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Abstract. 

Spectroscopy with HST/STIS Echelle and VLT/UVES of the bright QSO HE 0515-4414 (z cm = 1.73, B = 15.0) 
offers for the first time the opportunity to study the Lyman a forest in the redshift range 0.9 < z < 1.7 at a 
resolution < 10 kms -1 . The number density evolution of the Lyman a lines is well described by the power law 
approach dn/dz oc (1 + z) 7 . We derive 7 = 2.23±1.21 for the strong lines (13.64 < logA^Hi < 16.00), in agreement 
with the Lyman a forest evolution for z > 1.7. The expected slow-down in their evolution does not appear earlier 
than z ~ 1. For the weak lines (13.10 < log Am < 14.00) we find that the HE 0515-4414 data for 2 > 1 follow 
the trend with 7 = 1.40 known from z > 1.7 observations, i.e. we confirm the difference in evolution between 
weak and strong lines. We use the two-point velocity correlation function (TPCF) to search for clustering of the 
Lyman a lines, yet we detect no excess in the TPCF on scales up to 10 000 kms -1 . 
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1. Introduction 

The evolution of the Lyman a absorption lines per unit 
redshift dn/dz has been a subject of observational studies 
since many years. At high redshift [z > 1.6) the evolu- 
tion is steep with dn/dz oc (1 + z) 7 and 7 = 2 — 3 for 
strong Lyman a lines (i.e. log A'hi > 14; Rauch 1998[ ). 
Observations of the Lyman a forest at low reds hift with 
HST/FOS, at first in 3C273 flBahcall ct al. 199l|) , showed 
that the number of Lyman a lines was far in excess of the 
expected number according to an extrapolation from high 
z. As a fin al result from the HST QSO absorption line 
key project Weymann et al. (1998) found 7 = 0.1 — 0.3 
for z < 1.5. They also claimed a break in the evolution- 
ary behaviour from a steep evolutionary law for z > 1.6 
to a flat one for z < 1.5. The apparently rather abrupt 
break in the evolutionary law just at the transition from 
high-resolution optical data to low-resolution UV data im- 
mediately raised the suspicion that this behaviour might 



not be real. Owing to the insufficient FOS resolution 
(230 — 270 kms -1 ) the numbe r counts might be un der- 
estimated due to line blending ( Weymann et al. 1998 ). It 
is obvious, that this open question can be addressed only 
by UV observations of very bright QSOs at Echelle reso- 
lution. A further unsettled question which requires high- 
resolution UV spectra is the different evolutionary be- 
haviour of strong compared to weak Lyman a lines. It 
has been found from high- resolution optical spectra for 
z > 1.6 that apparently weak lines evolve much slower 



than strong lines (Kim et al. 1997). Is that also true for 
z < 1.5? 

In this paper we use the first UV spectra of the bright 
inter mediate redshift QSO HE 0515-4414 (z = 1.73, B = 
15.0, [Reimers et al. 1998j) at ST/S/Echelle resolution (10 
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km s ) for addressing the above discussed open questions. 



2. Observations 

HE 0515— 4414 was observed with STIS between January 
31 and February 2, 2000 with the medium-resolution NUV 
echelle mode (E230M) and a 0.2 x 0.2 aperture. The over- 
all exposure time was 31 500 s resulting in a typical signal- 
to-noise ratio of 10 depending on the order and on the po- 
sition within the orders. The resolution of the spectra is 
FWHM~ 10 kms -1 . The data reduction was performed 
by the HST pipeline completed by an additional inter- 
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order background correction and by coadding the separate 
subexposures. 

The optical spectra were obtained between October 
7, 2000 and January 3, 2001 using the UV-Visual 
Echelle Spectrograph (UVES) at the FLT/Kueyen tele- 
scope. The overall exposure time was 31 500 s. The slit 
width was 0.8 arcsec resulting in a spectral resolution 
of FWHM ~ 6 kms- 1 . After reduction by the UVES 
pipeline and conversion to vacuum baryocentric wave- 
lengths, the individual spectra were coadded and exhibit 
a S/N ~ 10 — 50 in the investigated spectral region. 

3. Data analysis 

The combined HST and VLT data provide the spec- 
tral range of the Lyman a forest of HE 0515— 4414 from 
z = 0.87 up to z — 1.73, the quasar's Lyman a emis- 
sion redshift. To avoid the proximity effect we exclude 
a region of about 5000 kms^ 1 from the quasar leading 
to an investigated spectral range A = 2278 — 3260 A or 
z = 0.87 — 1.68. We normalize the spectrum fitting poly- 
nomials to line-free regions and dividing the flux by this 
background continuum. 

The spectrum of HE 0515— 4414 is strewn with metal 
lines and lines of molecular hydrogen (H2) from its 
damped Lyman a system which has been studied in de- 



tail by de la Varga et al. (2000). The main difficulty in the 
line identification was the extraction of the Lyman a lines 
from this real H2 forest dominating the short-wavelength 
region of the HST spectrum and thus suggesting to ignore 
this region in our analysis of the weak Lyman a lines (see 
below). The molecular hydrogen of the DLA will be the 
topic of a forthcoming paper. 

In a first approximation, we detected about 400 lines 
as Lyman a candidates in the whole spectrum. These lines 
were fitted with the FITLYMAN code in the MIDAS pack- 



age (Fontana & Ballester 1995) using Voigt profiles con- 
volved with the instrumental profile. FITLYMAN adjusts 
three independent parameters per line by x 2 minimiza- 
tion. These parameters are the redshift of an absorption 
line z, its H I column density A/hi, and its Doppler param- 
eter b, comprising the thermal and the turbulent broaden- 
ing of the lines. The general fitting strategy, especially for 
blends, was to start with a single line and to add a further 
component to the ensemble if the \ 2 decreased with this 
second line. For each of the Lyman a candidates we calcu- 
lated the significance level S = — , where W denotes the 
observed equivalent width of the line and o~w the 1 a error 
of W implying both the fit error and the continuum error. 
With the selection criteria S ^ 1 and b )> 10 kms -1 , we 
reduced the original sample to 235 Lyman a lines. The 
full fit parameter list for all recognized hydrogen absorp- 
tion lines is available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr. 



4. Discussion 
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Fig. 1. Column densities of individually treated Lyman a 
lines versus column densities of simultaneous Ly a/Ly (3 
fits. We present the results for all detected pairs. 



4.1. Simultaneous fits with Lyman f3 

An important fitting constraint is given for the stronger 
Lyman a lines for which we tried to detect the higher or- 
der lines of the Lyman series followed by a simultaneous 
fit, if possible. We compared the Hi parameters inferred 
from Ly a profile fitting with the parameters determined 
with a simultaneous fit with the accompanying Ly [3 coun- 
terparts. Fig. 1 shows a rather small systematic effect in 
the column densities. We derive ( N^™ / N^ slc ) = 1.145 
for the mean proportion of the Ly a/Ly [3 fit column 
densities and the column densities based solely on Ly 
a fitting. The same holds for b (( 6 s i m / 6 s in g ic) = 1-03). 
This is in contrast to the recently suggested strong de- 
pendence of the fit parameters on the fit strategy (single- 
component fit versus simultaneous Ly a-Ly [3 fit; e.g., 
Hurwitz et al. 199$ |Shull et al. 2000|). Our 15 per cent 



deviation of the inferred column densities is at variance 
with the above cited studies. Obviously, the improved 
data quality leads to more consistent fit results. The re- 
maining discrepancy may be attributed to the unphys- 
ical assumption of Voigt profile fitting. Indeed, the in- 
terpretation of the line-broadening velocity as unresolved 
stochastic motions is probably an oversimplification (e.g., 
Levshakov fc Kegel 199?| ). 



4.2. Column density distribution function 

The differential column density distribution function 
/(A^hi) is defined as the number of Lyman a absorption 
lines per uni t column de nsity and per unit absorption dis- 
tance path ([Tytler 198%: 



/(JVhi) 



ANY,, AX, 



(1) 



For qn = 0, X(z) = \ [(1 + z) 2 - 1], so that Al = 
(1 + z)Az. For the spectral range investigated we get 
Az = 0.808, z = z = 1.278, and hence AX = 1.841. 



E. Janknecht et al.: First high-resolution study of the Lyman a Forest at 0.9 < z < 1.7 



3 



-10 
-1 1 

-12 



cn -13 

o 



-14 

-15 
-16 




13.00 < log N„ < 15.70: f(N„) « N„ 



1 3 



14 

log N„, 



16 



1 60 

140 

3 120 
c 

X> 

100 



60 



dato points 
best fit 

2 o confidence limits 



dn/dz = (111.60 ± 60.13) (1+z) - ' 1 a64 



1.4 



Fig. 2. The column density distribution function for the 

Lyman a lines. For the best fit only lines with log A Hl > Fig 3 The number density evolution of the Lyman a 



13.00 were included due to the limited S/N of the spectra. 



forest for weak lines (13.10 < N m < 14.00). The data 
points arc binned with Az = 0.1. The best fit was obtained 
by x 2 minimization. The dotted curves represent the 95 
% confidence band. 



4.3. Number density evolution 



Usually, /(Am) is fitted by a power law of the form 
/(JV H i) = AJV^f. The distribution function for all 232 
lines with log Nm < 15.70 is plotted in Fig. 2. The squares 
show the observed log / values, while the solid line rep- 
resents the best fit for all lines with log Ahi > 13.00. We 
have chosen this lower boundary because the distribution 
follows the power law down to this value. Assuming the 
general validity of the power law the sample is obviously 
not complete below log Am ^> 13.00. 

The best fit yields log A = 9.62 ± 0.58 and f3 = 1.61 ± 
0.04. Considering exclusively the STIS lines in the same 
column density range, the result is very similar: log A = 
9.47±0.59 and (3 = 1.60±0.04. In contrast, the slope of the 
distribution function for the UVES lines is much flatter 
(log A = 4.60 ± 0.72, = 1.24 ± 0.05). The distribution 
of the UVES lines can be approximated very well with a 
single power law for an even wider column density range 
(12.50 < log Am < 15.70) due to the better resolution of 
UVES. The flatter slope of the higher redshifted UVES 
lines indicates that stronger absorbers have evolved away 
faster than weaker ones. This will be discussed in more 
detail in Section 4.3. 

Our result is in accordance with other anal- 
yses in comparable redshift ranges. For example, 
Dobrzycki et al. (2002)| found /3 < 1.6 - 1.7, deriving the 
exponent from a curve of growth analysis. Hu et al. (1995) 
obtained (3 = 1.46 for 12.3 < log Am < 14.5, while 



Kim et al. (2001) determined (3 for various column den- 
sity ranges to 1.70 — 1.74. 



The evolution of the number density per unit redshift of 
Lyman a clouds can be well approximated by the power 
law 



An 

Az 



An 
d7 



(1 + ^) 7 , 



(2) 



where (gj) is the local comoving number density and the 
exponent 7 includes the cosmological evolution as well as 
the intrinsic evolution of the absorbers. For a non-evolving 
population in the standard Friedmann universe with the 
cosmological constant A = and with qo = the exponent 
becomes 7 = 1. 

Because there exist a lot of hints that weaker and 
stronger Lyman a clouds evolve differently, it is conve- 
nient to distinguish between two subsamples. We use the 
column density range 13.10 < log Am < 14.00 for the 
weak lines, while the interval 13.64 < log Am < 16.00 de- 
fines the strong ones. These boundaries are frequently used 



in the literature ( 


Weymann et al. 1998; 


Kim et al. 2001 


Dobrzycki et al. 2002 


) allowing a direct comparison. 



In Fig. 3 and Fig. 4 we present the line numbers per 
unit redshift plotted over the redshift for the weak and 
for the strong lines, respectively. While for the strong 
lines we could exploit the whole wavelength region, we 
omit the spectral range A ^ 2415 A for the weak lines 
to avoid misidentification due to the H 2 lines, retaining 
an effective redshift range 1.0 ^ z < 1.7. The diagrams 
show the data points and the best fit for which we obtain 
7 = 0.01 ± 0.64 for 13.10 < log A Hl < 14.00, suggesting 
that there is little evolution in the weak lines in the red- 
shift interval 1.0 < z < 1.7. Considering a broader redshift 
range (Fig. 5, upper part) demonstrates that our STIS 
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Fig. 4. Same as in Fig. 3, for 13.64 < N Hl < 16.00. 



data follow the earlier optical observations (z > 1.6), i.e. 
the number density evolution is consistent with 7 = 1.40 
over the whole redshift range 1 < z < 3.7. It should be 
noted that 7 = 1.40 lies within the 2er confidence band 
of our data points (see Fig. 3). The transition to a flat 
(7 = 0) evolution curve probably occurs around 2=1. 

In contrast, the strong Lyman a lines show a steeper 
gradient in the evolution diagram (Fig. 4). We detect an 
obvious correlation between the evolution and the line 
strength, i.e., the high column density absorbers evolve 
with 7 = 2.23 ±1.21. This disagrees with the results 



of [Pcnton ct al. (2000)| and pobrzycki ct al. (2002)| who 
found no or only marginal evidence for a different evo- 
lution, respectively. 

We cannot recognize a slow-down in the evolution of 
the stronger absorbers. Therefore, we conclude that this 
break does not occur earlier than at z ~ 1 rather than at 



z ~ 1.5 — 1.7 as pre v iously claimed (Impey et al. 1996 



Weymann ct al. 1998j Pobrzycki ct al. 2002[ ). The large 



spread of our data points result from the poor statistics of 
a single line of sight. For example, omitting the two out- 
liers in the plot for the strong lines lying beyond the 95 
% confidence limit our result (7 = 1.88 ± 0.57) becomes 
more robust. 

The lower panel of Fig. 5 demonstrates the difference 



between the Weymann et al. (1998) data points and our 
values. Indeed, the former ones do not indicate any change 
in the evolution until z ~ 1.5, while the results of HE 
0515—4414 suggest a change of the slope at much lower z. 

4.4. Two-point velocity correlation function 

To study the clustering properties of the Lyman a forest, 
we introduce the two-point velocity correlation function 
£ (Av) where the number of observed line pairs in a given 
velocity separation bin Av, n Q bs, is compared with the 
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Fig. 5. The number density evolution of the Lyman a for- 
est, comparison of different studies. Shown are the results 
for the column density range iVm = 10 1310_14 00 cm~ 2 
(upper panel) and JVhi = l0 13 - 64 ~ 16 00 cm~ 2 (lower panel). 
The filled circles are from our study, for the other sym- 
bols see legend. It should be noted that the data points 
from |Pcnton et al. (2000)| and |Savaglio et al. (1999)| in the 



lower panel apply to iVm > 10 cm 



number of expected pairs n C xp determined in the same 
velocity difference bin in a randomly produced spectrum: 



e(A») 



n obs (Av) 
(Av) 



(3) 



^cxp 



Having two absorbers at z\ and zi the velocity splitting 
in the rest frame at mean redshift reads 



Av = 



c(z 2 - Zl) 



zi+z 2 
2 



(4) 



We derived n oxp from Monte Carlo simulations distribut- 
ing our line sample - in accordance with the determined 
number density evolution - randomly over the spectrum, 
counting the line pairs for varying velocity splittings and 
repeating this procedure 1000 times. 

Various studies have found at least marginal evidence 
for clustering on scales Av < 500 kms - . For exam- 
ple, |Cristiani et al. (1995)| derived £ = 0.89 ± 0.18 and 



£ = 1.02±0.26 (for two different QSOs) for log N Kl > 13.8 
and a redshift interval < z < 3.66 at a velocity separa- 
tion of Av = 100 kms -1 , while [Kim ct al. (2001)| found 



£ = 0.4 ± 0.1 for logiVHi > 12.7 and 1.5 < z < 2.4 
at the same velocity splitting. Pcnton ct al. (2000)| and 
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Fig. 6. Two-point correlation function £ (Aw) (solid line), 
in 100 kms -1 bins, for higher column density absorbers 
(13.64 < logiVm < 16.00). Dashed and dot-dashed lines 
represent the la and 2a Poisson errors, respectively. 



Kulkarni et al. (1996) also detected a weak clustering sig- 
nal in their data. In addition, the occurence of clus- 
tering on different scales is also supported by a num- 
ber of theoretical studies (e.g., |Cirkovic fc Lanzetta 200C 
Pando et al. 19981). 



We have examined the clustering separately for the 
weak and strong lines (the latter is shown in Fig. 6). Up 
to a velocity separation of Av = 1000 kms -1 we detect 
no signal exceeding the 2a level and only a marginal one 
above the la level. For example, for the strong lines we 
derive £ (Av — 50kms _1 ) = 1.14 with a la poisson er- 
ror of 1.13 in accordance with Fig. 6 . For log iVm > 13.8 
we obtain £(50kms _1 ) = 1.63 (la poisson error 1.60) 
and £(150kms _1 ) = 1.23 (la poisson erro r 1.36), so 
we cannot confirm the |Cristiani et al. (1995) result. For 
logiV H i > 12.7, £(50 kms -1 ) = 0.10 (la poisson error 
0.16) and £(150 kms" 1 ) = -0.01 (la poisson error 0.15) 
not being consistent with the weak clustering seen by 
Kim et al. (2001) in a comparable redshift range. There 



is also no excess in £ on much larger velocity separations 
(up to 10000 kms -1 ) in our data. We do not find any 
dependence of £ on the redshift either. Thus we cannot 
confirm the result of weak clustering of the Lyman a lines 
as suggested in the above cited studies. However, the small 
number of lines in a single line of sight does not allow a 
final assessment. 

5. Conclusions 

We have analyzed the evolution of the Lyman a forest 
in the redshift range 0.9 < z < 1.7 using combined high- 
resolution HST/STIS and VLT/UVES data. The main re- 
sults are summarized as follows: 



1. The recently claimed dependence of the hydrogen fit 
parameters on the fit strategy (single Ly a fit or si- 
multaneous Ly a-Ly [3 fit) cannot be confirmed by our 
results. 

2. The column density distribution of our Lyman a line 
sample can be approximated over several orders of 
magnitude (13.00 < log TVjji < 15.70) by a power law. 
The (negative) slope is described by (3 — 1.61 ± 0.04 
for ~z ~ 1.3, consistent with recent results. 

3. The evolution of strong and weak lines is distinctly dif- 
ferent. The high column density (13.64 < logA^i < 
16.00) absorbers evolve according to (1 + z) 7 with 
7 = 2.23 ± 1.21 for 0.9 < z < 1.7, and the expected 
slow-down in the evolution does not appear down to 
z ~ 1 . The evolution of the weaker lines over the same 
redshift range is consistent with 7 = 1.40, thus we 
have a continuation of the trend seen at higher red- 
shifts. Again, the transition to non-evolution probably 
occurs around z = 1. More lines of sight are necessary 
to confirm our results. 

4. We detect no significant clustering of neither the weak 
nor the strong Lyman a lines on scales up to 10000 
kms -1 , so we are unable to confirm most of the pre- 
vious studies reporting a weak clustering signal for 
Av < 150 kms" 1 . 
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